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£\\ [ Abstract 

It is shown that hypothetical neutrino-majoron coupling can suppress neutrino flavor 
oscillations in the early universe, in contrast to the usual weak interaction case. This 
\q I reopens a window for a noticeable cosmological lepton asymmetry which is forbidden 

\Q . for the large mixing angle solution in the case of standard interactions of neutrinos. 

o 

O 

1 1 Introduction 

Cosmological lepton asymmetry is not directly measurable, in contrast to baryon asym- 
metry, but may be observed or restricted through its impact on big bang nucleosynthesis 
(BBN), large scale structure formation, and the angular spectrum of the cosmic microwave 
• background radiation (CMBR), for a review see e.g. Ref. [1]. At the present time the best 

bounds follow from the consideration of BBN. Primordial production of light elements is 
especially sensitive to the value of asymmetry between electronic neutrinos and antineu- 
trinos since they directly influence the neutron-proton transformations in weak reactions 
v e n <-> e p and v e p «-> e + n. The bound on the chemical potential of electronic neutrinos 
obtained in Ref. [2] reads |£ e | = \fi e /T\ < 0.1. The bound on muonic or tauonic asymmetry 
is noticeably weaker because z/« or v T produce an effect on BBN only through their impact 
on the cooling rate of the universe and the degeneracy of these neutrino flavors is equivalent 
to an addition of 



X 

5-H 



7T / V 7T 



(1) 



massless neutrino species at the BBN epoch. If the latter is bounded by AN U < 1 (safe 
bound) then |^ jT | < 1.5. For less conservative bound, AN V < 0.4 [3], the dimensionless 
chemical potentials should be below 0.9. Still, asymmetry of the order of unity would have 
noticeable cosmological effects on large scale structure formation and CMBR. The limits 
quoted above are valid if the chemical potential of only one kind of neutrino is non-zero. If a 
conspiracy between different chemical potentials is allowed, such that a positive effect of 
or £ r is compensated by £ e or vice versa, the limits would be somewhat weaker. According 
to Ref. [4] they are: |£ e | < 0.2 and |£ M)T | < 2.6. 1 

1 We are not concerned in this paper with the origin of such large lepton asymmetry. Thus far, many 
cosmological scenarios that explain both the small observed baryon asymmetry and a possible large lepton 
asymmetry simultaneously, have been proposed [5]. 
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The bounds on chemical potentials of v„ and v T can be significantly improved because 
of the strong mixing between different neutrino flavors [6]. This mixing gives rise to the 
fast transformation between v e , v^, and v T in the early universe and leads to equilibration 
of asymmetries of all neutrino species. Thus the BBN bound on any chemical potential 
becomes essentially that obtained for v e [7] (see also the papers [8]): 

\^,r\ < 0.07. (2) 

In this case the cosmological impact of neutrino degeneracy would be negligible. However, a 
compensation of £ e by other types of radiation is still possible, which was studied in Ref. [9]. 
The authors of Ref. [9] put the constraints on £ e and AiVj, based on WMAP data and BBN: 

-0.1<£ e <0.3, -2<AJV„<5. (3) 

It is interesting to see if one could reasonably modify the standard model to allow large 
muonic and/or tauonic charge asymmetries, together with a small electronic asymmetry, 
to avoid conflict with BBN. This is the aim of this work. A natural generalization is to 
introduce an additional interaction of neutrinos with massless or light (pseudo)Nambu- 
Goldstone boson, majoron [10]. The idea to invoke majoron to modify neutrino oscillations 
in primeval plasma was discussed in Refs. [11, 12]. In these papers, two different mechanisms 
which could block or suppress oscillations between active and hypothetical sterile neutrinos 
have been considered. We have, however, some concerns about the validity of their results 
which we will discuss in the next section. Let us note that in this paper we consider an 
impact of neutrino majoron interactions on the oscillations between active neutrinos and 
not on active-sterile oscillations, as is done in the mentioned above papers. 



2 General discussion and approximate estimates 

Neutrino oscillations may be suppressed in medium if the interaction of neutrinos with the 
medium is sufficiently strong or in other words the refraction index, n (or effective potential, 
V = E(n — 1)) of neutrinos is very large. Correspondingly the mixing angle in medium 6 m 
becomes negligible: 

s 2 

tan 2 29 m = (4) 

(4 + V/5E) 2 + r 2 /(4<5£ 2 ) 

where S2 = sin 26, C2 = cos 26, 6 is the vacuum mixing angle, 5E = 5m 2 /2E, E is the 
neutrino energy, and T is the rate of neutrino interaction with medium given by 

_ 80(1 +gl+gl)a 2 F ET* 

where T is the temperature of the primeval plasma, Gf = 1-166 • 10 -11 MeV~ 2 is the Fermi 
coupling constant, gi = sin 2 6\v ±1/2 and gn = sin 2 6\v where sin 2 8\y = 0.23 and the sign 
"+" stands for u e and "— " stands for u^ T . The derivation of these equations can be found 
e.g. in lectures [13]. 

The diagonal components of the effective potential, created by the standard weak inter- 
action, for the active neutrino species are given by [14]: 

= ±C lV ^G F T 3 - C^^-^ (6) 
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where a = e,fi,T labels the neutrino flavors, a = 1/137 is the fine structure constant, and 
the signs "±" refer to neutrinos and antineutrinos respectively. The first term arises due to 
a possible charge asymmetry of the primeval plasma, while the second one comes from the 
non-locality of weak interactions associated with the exchange of W or Z bosons. According 
to Ref. [14] the coefficients Cj are: d « 0.95, Cf « 0.62 and C£ T « 0.17 (for T < m M ). 
These values are true in the limit of thermal equilibrium, otherwise these coefficients are 
some integrals from the distribution functions over momenta. The charge asymmetry of 
plasma is described by the coefficients r/ a ) which are equal to 

ri (e) = 1f)v e + + Vvr +rie~ W 2 ( for v e) , (7) 
r]V) = 2^ + n„ e + r] UT - rj n /2 (for , (8) 

and rf T ^ for v T is obtained from Eq. (8) by the interchange n <-> r. The individual charge 
asymmetries, rjx, are defined as the ratio of the difference between particle-antiparticle 
number densities to the number density of photons with the account of the 11/4-factor 
emerging from the e + e~-annihilation: 

Vx = (N x -Nx)/N 1 . (9) 

For "normal" values of charge asymmetry, i.e. \rj\ ~ 10 -9 the charge asymmetric term 
in the potential (6) is subdominant but if the asymmetry is of the order of unity, then the 
ratio V aa /5E becomes huge: 

V^/SE « 2 • 10- 11 r]ET 3 (5m 2 )- 1 (MeV) -2 « 20r/ (T/MeV) 4 (eV 2 /<5m 2 ) (10) 

and e.g. for 5m 2 = 10 -2 eV 2 and rj ~ 1 the mixing angle in matter would be suppressed 
more than by three orders of magnitude in the BBN range of temperatures, T ~MeV. 

However, this result is valid for mixing between active and sterile neutrinos and is not 
true for active-active mixing. In the last case the effective potential has large off-diagonal 
matrix elements [15], V a b, with a / b, which compensate the suppression induced by the 
diagonal terms. This is the reason for non-suppressed oscillations between active flavors 
and for the equilibration of all leptonic asymmetries [7]. 

Possible additional interactions of neutrinos with majoron can be described by the La- 
grangian: 

tint ~ X{9ab^Cu b + h.c) (11) 

where x is the operator of the majoron field, C is the matrix of charge conjugation, and 
g a b are the coupling constants. We will assume for simplicity, though it is not necessary, 
that only the flavor-diagonal coupling is non- vanishing, g ab ~ g a a^ab- We will consider the 
coupling matrix with more general form later. Let us also assume that the majoron is 
so light that its decay and inverse decay are not essential at the BBN epoch. Neutrino- 
(anti) neutrino scattering through majoron exchange gives rise to a contribution to neutrino 
effective potential (proportional to forward scattering amplitude) which can be estimated 
as 

~ g aa g b bT. (12) 

The constants g aa should satisfy several constraints to make the mechanism operative. 
Firstly, the diagonal part of the potential Va$ should be larger than the weak potential V^*") 
given by Eq. (6) , while its off-diagonal components must be much smaller than the diagonal 
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ones, « Vaa , that is, the flavor symmetry in the neutrino-majoron interactions 

should be strongly broken. These two conditions would ensure suppression of flavor changing 
oscillations between active neutrinos. To realize these conditions the following constraints 
should be imposed: 



Here a labels |iorr and the coupling of majoron to u e is assumed to be much weaker than 
those to Vn and/or v T \ T ae is the temperature at which u a are effectively transformed into v e 
in the standard theory. According to calculations of Ref. [7] it takes place around T ae = 1 
MeV for the large mixing angle solution to solar neutrino deficit, while the transformation 
between and v T in the standard theory takes place somewhat above 10 MeV. In fact, as 
we see in what follows from numerical calculations, the more accurate bound is much less 
restrictive than this simple estimate (see Eq. (49) or Figs. 6 and 7 below). 

The second inequality (13) implies that the off-diagonal components of the effective 
potential are suppressed in comparison with the diagonal ones and the oscillations remain 
blocked. This is not so in the case of the standard weak interactions which are flavor 
symmetric and the large value of the denominator due to large V aa in Eq. (4) is compensated 
by the same large off-diagonal components of the effective potential. 

Secondly, the coupling constants g aa should not be too large, otherwise flavor non- 
conserving reactions of the type v e v a <-> T> e u a (or similar) would lead to equilibration of all 
leptonic charges. To avoid that the rate of these reactions, T ea ~ a ea T 3 , should be smaller 
than the cosmological expansion rate H ~ T 2 /mpi, where mpi = 1.221 • 10 22 MeV is the 
Planck mass. Thus, to suppress e — \x or e — r transformation through direct reactions one 
needs 



This conditions should be satisfied for temperatures above the BBN range, i.e. T > 1 MeV. 
Similarly, if we require that u a v a <-> v a v a should not occur efficiently, the coupling constants 
must satisfy a similar inequality with g ee replaced with g aa . 

There are quite strong limits on possible coupling of majoron to neutrinos which follow 
from astrophysics [16, 17]; discussion and references to earlier works can be found e.g. in 
the book [18]. Astrophysics allows either very small or quite large coupling constants. The 
former is quite evident, while the latter appears because strongly interacting majorons, 
though efficiently produced inside a star, cannot propagate out and carry away the energy, 
thus opening a window for large values of the coupling. It is not so for the coupling 
to v e because the latter is bounded from above by the data on double beta decay [19], 
<7 ee < 3 ■ 10~ 5 . Together with the supernova bounds, the upper limit is shifted down to 
g ee < 4 ■ 10~ 7 [17], with a small window around (2 — 3) • 10~ 5 . So we assume in the 
following that g ee < 4 • 10~ 7 . For fi or r the allowed regions are: g aa < (3 — 5) • 10 -6 or 
9aa > (3 — 5) • 10 -5 . Evidently the conditions specified above can be satisfied. As reference 
values we may take g ee = 1CT 7 and g aa = 5 • 10~ 6 which satisfy all constraints presented 
above and would lead to a suppression of the transformation of or v T into v e , thus 
permitting a large muonic or tauonic asymmetry combined with a small electronic one at 
BBN. Note that it is not necessary to assume a large value of g aa ^ 10~ 5 to accomplish this. 
As shown later, such a large coupling constant would change the scenario into something 
more complex. 

If majoron is lighter than the mass difference of neutrinos then a heavier neutrino would 
decay into a lighter one and majoron. This process might leave traces in the flavor ratios of 




(13) 




(14) 
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the high-energy neutrinos from distant astrophysical sources [20], which can be detected by 
e.g., IceCube [21]. This could be potentially sensitive to very small values of the neutrino- 
majoron coupling. Existing bound [22] on life-time/mass of decaying neutrinos based on 
the solar neutrino data, r/m > 10~ 4 sec/eV, is too weak to be essential for the mechanism 
discussed in the present paper. 

In the subsequent sections we will present more accurate calculations but before turning 
to a closer examination of the scenario, let us make a few comments on earlier papers where 
majoron suppression of neutrino oscillations in the early universe have been considered. In 
Ref. [11] it was assumed that there exists a coupling of an active and sterile neutrinos to 
majoron. Let us note that in the version of the theory that we consider here no sterile 
neutrinos are introduced: the negative helicity states are identified with neutrinos, while 
the positive helicity states are identified with antineutrinos. The coupling considered in 
Ref. [11] has the form: 

£as = Was XWs + h.C. (15) 

where u a is an active neutrino flavor and v s is a sterile one. The authors argued that the 
effective potential induced by the interactions of active and sterile neutrinos with majoron 
strongly suppressed the oscillations. However they omitted off-diagonal V^ s -term in the 
effective potential which might invalidate their result. Possibly this mechanism of oscillation 
suppression may operate in a more complicated version of the model. 

In Ref. [12] a different mechanism of neutrino oscillation blocking by majoron has been 
suggested. According to the equation of motion, the total leptonic current including majoron 
and neutrino contributions is conserved: 

fD 2 X + = (16) 

where is the covariant derivative in cosmological background, / is the vacuum expec- 
tation value of the Higgs field responsible for the spontaneous breaking of leptonic ?7l(1)- 
symmetry and J M is the leptonic currents of fermions. In spatially homogeneous case the 
solution to this equation is 

, 1 3 



nil) v Xin -•/?„//) +•/*(*)//, (17) 

where a(t) is the cosmological scale factor and the sub-index "in" indicates initial values. We 
assume that Xin 

= 0. Hence the authors of Ref. [12] concluded that x = VB^/f, where r) B = 
6- 10~ 10 is the baryon asymmetry of the universe. The contribution to the neutrino effective 
potential from this term is 5V X = x/ f and, according to [12], with f 2 = (5 — 9) GeV 2 the 
potential would be strong enough to suppress the oscillations. However, according to the 
usual scenarios, baryon asymmetry could be created at much larger temperatures, much 
higher than this scale. Thus, after spontaneous symmetry breaking, leading to creation 
of majoron, J* remains constant in comoving volume and x = 0. A non-zero x might be 
created if lepton asymmetry was generated at spontaneous breaking of leptonic Ul(1) or 
by the neutrino oscillations themselves but both these cases need further and more detailed 
investigation. So in our considerations we assume that the mechanism of Ref. [12] is not 
operative. 



3 Effective potential induced by neutrino-majoron interac- 
tions 

In this section we derive the effective potential induced by the neutrino-majoron interac- 
tions, which is one of essential ingredients of our paper. The relevant part of the Lagrangian 
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is given by 



C 

C>int 



a 

\x(gab ® T a <J2®b + 9* ab > 
\x (gab v T a Cv h + g* ab 4 C <) > 



(18) 



(19) 



where $ a and v a are two-component and four-component representations of neutrino of 
flavor o, respectively. They are related to each other as = (<3?^,0) T in the chiral rep- 
resentation (see Appendix A for notations). Here and hereafter u a is taken to be the 
left-handed field. In the following we neglect the effect of small masses of neutrinos, and 
treat them as massless fields. In this limit, there are three diagrams 2 that contribute to 
the effective potential (see Figs. 1 and 2 ). 

Let us first calculate the effective Hamiltonian corresponding to neutrino- (anti) neutrino 
scattering processes. For that purpose we integrate out the majoron field using its equation 
of motion, and substitute the solution into interaction Lagrangian (19). Then we quantize 
neutrino fields perturbatively. The details of the derivation can be found in the Appendix A. 
The effective Hamiltonian, which describes the neutrino- (anti) neutrino scattering processes, 
is given by 
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/ 



dp <iq dr ds 



5^ (p + q + r + s) e- iEtott F(p, q)F(r, s) 
epegEpE^l - e p e q cos 6> pq ) 



(20) 



with 



E P = |p| , 
Etot = CpEp + tqEq + e r E r + e s E s , 

COSflpq = i i i i - 

IpI |q| 

^(P.q) = 9abVa{p) T Cv b (q) +SW ) ffe(-q) t O*(-p), 
where dp = d 3 p/(27r) 3 . The momentum expansion of free neutrino field is 



n{x) = J dp ^( P ) e ^°* +i P x , 

= J dp (a t (p)u p + &i(-p)t v _ p ) e~^ t+ ^ 



(21) 



(22) 



with p° = e p \p\ = ± |p|. The sign of e p is chosen to reproduce positive or negative energy 
solution. Here aj(p) (^(p)) is the annihilation (creation) operator for negative (positive)- 
helicity neutrinos with momentum p, u p (v p ) represents left-handed Dirac spinor of negative 
(positive)-helicity neutrinos. With thus obtained effective Hamiltonian, the effective poten- 
tial can be calculated by the technique of Sigl and Raffelt [23]. Although the calculations 
are lengthy, the procedure is straightforward. The contribution to the effective potential 
for neutrinos with momentum p is given by 



ab 



(iq 



4 IpI |q| 



(Pq + Pq) 9 



ab 



(23) 



Note that there are other diagrams if small majorana masses are taken into account. However, the 
amplitudes of such diagrams are suppressed by powers of m v /E u <JC 1 for rclativistic neutrinos. 
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Figure 1: s- and t-channel diagrams of neutrino- (anti)neutrino scattering through majoron 
exchange. Both diagrams contribute to the off-diagonal component of the effective potential. 



where the density matrices are defined as 

[a\(p) a z (p')) = (2vr) 3 ^)(p-p')[p p ; 



4(p)Mp': 



(2^) 3 5 (3) (p 



P 



(24) 



If we take the coupling constant matrix g a b to be diagonal, the effective potential becomes 



v p 



ab 



9aa9bb / dq 



1 



4|p| |q| 



(Pq+Pq 



ab 



(no summation), 



9aa9bb] 



This result confirms the rough estimate of the effective potential in the previous section. 

The neutrino-majoron scattering shown in Fig. 2 also contributes to the effective poten- 
tial and can be calculated similarly. The result is 



yiyx) 



ab 



r/q 



/x(q) 

4|p| |q| 



9 ] 9 



(25) 



where /v(p) is the number density of majorons with momentum p. Note that this expression 
vanishes if the abundance of majorons is negligible in comparison with their equilibrium 
value, i.e. f x (p) ^ /x (p)- Thus the complete effective potential for neutrinos with 
momentum p induced by interactions with majorons is given by 



y(x) 
v p 



ab 



dq 



1 



4|p| |q| 



(26) 



where 1 is the unit matrix in the flavor basis. 



4 Possible role of neutrino-majoron reactions. 

As we have already mentioned reactions between neutrinos with an exchange of majorons 
do not conserve leptonic charge and if they are efficient, lepton asymmetry of the universe 
would be completely destroyed before BBN. One of the dangerous processes is 

2u <-> 29. (27) 

We assume for simplicity that the diagonal coupling constants \g a a\ are larger than non- 
diagonal ones and that one flavor coupling dominates, e.g. \g TT \ > > \g e e\- In this case 
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Figure 2: Diagram of neutrino- majoron scattering through neutrino exchange. If majorons 
are not abundant in the primeval plasma, this diagram gives negligible contribution to the 
effective potential. 

only three diagrams presented in Fig. 3 contribute into the reaction (27). The amplitude 
corresponding to any of these three diagram is equal to \g\ 2 and thus the complete amplitude 
squared is simply | A(2f «-> 2t>)| 2 = 9\g\ 4 (see Appendix B). Here and in what follows we 
omit the sub-index oo at the coupling constant g aa . With the known amplitude of the 
L-nonconserving reaction (L is the leptonic charge), the kinetic equation can be written as 

— xd x f v {yi) = - ■ / F(f) 5 { >{ Vl + 2/2 - 2/3 - 2/4 + I^J 28 

x z 2 2»ir b xyi J 2/2 2/3 2/4 

where / is the neutrino distribution function, (yi,yi) = (Ei/T,pi/T), x = (T/MeV) _1 , T 
is the temperature of the cosmic plasma, (iJ/MeV) = h/x 2 is the Hubble parameter with 
h ~ 4.5 • 10~ 22 , J^j] is the contribution to the collision integral from reactions that conserve 
leptonic charge, and 

F(f) = [l - / P (yi)][l - HyzWMMvA) - [l - fu(y 3 )][l - MviftMvdUfa). (29) 

We assumed that the temperature evolved according to T = —HT. The combinatorial 
factor 1/2 in front of the r.h.s. comes from two factors 1/2! and one factor 2 because two 
identical particle are produced (see e.g. Ref. [24]). 

The collision integral in Eq. (28) can be easily evaluated in the limit of Boltzmann 
statistics, when / C 1 and under assumption of kinetic equilibrium, so the distribution 
functions have the form / = exp(— y + £), where £ = n/T is the dimensionless chemical 

(el) 

potential. To avoid the contribution of I^j to the equation governing the evolution of £ 
we take the difference of Eq. (28) and analogous equation for f u and integrate it over (i 3 yi. 
The conti 
the form: 



The contribution of 1^9, to this integrated difference vanishes. The resulting equation takes 



where we have used the formula presented in Appendix C. It is also assumed above that 
£ = — £. Otherwise the contribution from L-conserving ^-annihilation would not vanish in 
the collision integral. If £ > 1 we may neglect exp(— £) and solve equation (30) as 

eXpK(x)] = l + 9M 4 zexpf (0)]/(128^) (31) 
The lepton asymmetry would not be destroyed by reaction (27) if 

\g\ < 2- 10- 5 T 1 / 4 exp[-£(0)/4]. (32) 
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Here (and in what follows) T is given in MeV. Quantum statistics (Fermi) effects somewhat 
weaken the bound. Their presence does not allow to solve kinetic equation analytically. 
Numerical estimates give a similar constraint as shown in Fig. 5. 

Muonic or tauonic charge asymmetry should be preserved until the annihilation v T ^v T ^ 
into e + e~-pairs is frozen. Otherwise we will return to the standard scenario with zero lepton 
asymmetry. According to the estimates of Ref. [13] made in Boltzmann approximation, the 
freezing temperature of the annihilation is Tj « 5.3 MeV. Fermi corrections would make 
its value slightly higher. If muonic and/or tauonic charge asymmetry were erased by the 
oscillations or L-nonconserving reactions with exchange of majorons below Tf then the total 
number density of v T plus P T (or + P^) would be conserved in the comoving volume and 
their distribution would be given by 

U a = fv a = \exv(ElT-t a )\- x (33) 

where a = /j,,t. So effectively the energy density of these neutrinos would be the same 
(up to a numerical factor of order of unity) as the energy density of the usual degenerate 
neutrinos. If the mixing of v T or with v e would not change the distribution of the latter 
then BBN would allow a large values of £ a and the degeneracy of v a may lead to noticeable 
cosmological effects. 

We make a simplifying assumption of absence of majorons in noticeable amount in the 
primeval plasma. If majorons would be in equilibrium we should take into account their 
(4/7)-contribution into the number of effective neutrino species at BBN and modification of 
neutrino effective potential due to forward elastic scattering v\ — ► vx- Even if majorons arc 
abundantly produced, the main conclusion of our paper remains unchanged but numerical 
results would be somewhat different. The production of majorons could proceed through 
the reaction v + P — > 2%. The Feynman diagrams describing this process are presented in 
Fig. 4. The amplitude squared of this process is 

,2 I |4 ( Pi ' Ps , Pi -P4 
\M- mv \ = \g\ 1 

\Pi - Pa Pi- Ps 

where Pim an d £>3,4 are four momenta for (anti) neutrino in the initial state and majorons 
in the final state, and they satisfy the momentum conservation condition, P1+P2 = P3 + Pi- 
Note that the collision integral with this amplitude squared involves an IR logarithmic di- 
vergence, therefore we need to input a lower cutoff scale. Taking into account the thermal 
effects, the dispersion relations for the majoron and the neutrinos change from those at 
zero temperature. Especially, they obtain finite thermal masses, which provide the desired 
infra-red cutoff. For \g\ ~ 0(1CP 5 ), the thermal effect of the neutrino-majoron interac- 
tion dominates over that due to the electroweak interaction, so the IR divergent part is 
regularized as ~ \og{AE 2 / m\) with tut ~ \d\T. 

Kinetic equation governing the production of majorons is similar to Eq. (28), and has 
the form: 

h | g | 4 log(4£ 2 /r4) exp(-y) 

27jt3x y +- W 

where we omitted elastic reactions. Integrating this equation we can find for the ratio of 
the number density of % to the equilibrium value: 

% _ M 4 log(4£ 2 /"4) x ^ 



n 



N) " 2 8 C(3)tt 3 /i 
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Figure 3: s-, t- and u-channel diagrams of 2v <-> 2v through majoron exchange. 
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Figure 4: t- and u-channel diagrams of v + v *-* 2\ through neutrino exchange. 



Demanding that n x < n x we obtain 

\g\ < 2.1 • lO^T 1 / 4 , (36) 

which gives a constra int similar to Eq. (32). However, if we take into account quantum 
statistics, the bound becomes slightly weaker for large charge asymmetry of neutrinos, as 
shown in Fig. 5. 

So, to summarize, the muonic or tauonic lepton asymmetry would be erased if the 
neutrino-majoron interaction is sufficiently strong, i.e. \g\ ^ 10 -5 . If we restrict ourselves to 
a scenario that the majoron-neutrino interactions suppress neutrino oscillations and thereby 
keep the large lepton asymmetries intact, the largest coupling constant should be smaller 
than ~ 10~ 5 . However this does not necessarily mean that the conspiracy between the 
speed-up effect of and a shift of (3 equilibrium due to £ e is impossible for \g\ ^ 1CP 5 . 
In fact, if the L-nonconserving interactions became efficient only after muonic and tauonic 
neutrinos decoupled from the thermal plasma, then the neutrino degeneracy would still be 
maintained in muonic or tauonic neutrinos resulting in their larger energy density at BBN. 
Majorons, as well, can contribute to additional relativistic degrees of freedom at BBN. We 
will discuss these issues below. 



5 Neutrino oscillations in primeval plasma 

As is established now, neutrino mass eigenstates are related to the flavor eigenstates through 
the orthogonal matrix (we neglect a possible CP violation): 

Va = U aj Vj (37) 

where a = e,fi,T and j = 1,2,3. We have chosen the parameters so that in the limit of 
small mixing v e ~ v\, Va ~ V2, v T ~ ^3- However, the mixing between neutrinos are known 
to be large and they cannot be taken as a dominant single mass eigenstate. 
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Figure 5: Upper bounds on \g\ obtained from the requirement that L-violating reactions, 
2u *-* 2u and v + u <-> 2%, are out of equilibrium. Numerical and analytical bounds for the 
former reaction are shown as solid and long-dashed lines, respectively. The short-dashed 
line represents the constraint obtained numerically for the latter reaction. 

Kinetic equations for the density matrix of oscillating neutrinos can be written as [25, 23]: 

ip=[nw,p}-i{n^,p} (38) 

where the first commutator term includes the vacuum Hamiltonian and the effective po- 
tential of neutrinos in medium calculated in the first order in Fermi coupling constant, Gp 
(weak interaction part) and in the second order in the coupling to majoron. The latter 
is given by Eq. (26) and does not contain off-diagonal terms if only one flavor coupling 
constant dominates. The weak interaction potential contains off-diagonal terms of the same 
order of magnitude as diagonal ones because of universal coupling of W and Z bosons to 
all neutrino flavors. Their explicit expressions can be found e.g. in Ref. [7]. In what follows 
we will skip upper indices (1) and (2) at TL. 

The second anti-commutator term in Eq. (38) describes breaking of coherence induced by 
neutrino scattering and annihilation as well as neutrino production by collisions in primeval 
plasma. It includes imaginary part of the Hamiltonian calculated in the second order in 
Gf and the fourth order terms in neutrino- majoron coupling g related to neutrino- neutrino 
scattering through majoron exchange. If majorons were abundant in the primeval plasma 
the processes of neutrino-majoron scattering should also be included. 

Since the "atmospheric" neutrino mass difference is much larger than the "solar" one, we 
may simplify the problem considering the process in two steps. First, oscillations between 
v T and Vn were switched-on. This started at temperatures about 10-15 MeV and might 
lead to equilibration of tauonic and muonic charge asymmetries if these asymmetries were 
not erased by reactions (27). Later at T < 3 MeV, oscillations would start between v e and 
ty , which is a certain mixture of cv T + sv^. This process is potentially dangerous for BBN 
because it could change number density and spectrum of v e and v e . 

If mixing is effective only between two neutrinos, the density matrix is 2 x 2 and kinetic 
equations for its components have the form: 

iHx8 xP ^ = H^Pr^-Hr^r-il^ (39) 
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iHxd xPTT = Hr^r-H^Pr^-il^ (40) 
iHxd x p^ T = (Hup - TCrr) Pixr + H^T (Ptt - Pwi) - i^iirPjXT (41) 

where H and x are defined below Eq. (28), TL ab = ^ab+Sj=2,3 U a jU b jm 2 /2E, U a j are matrix 
elements of the mixing matrix; in the case under consideration U^ 2 = U T % = cos 8 = c and 
[7^3 = — U T 2 = sin 8 = s. The coherence breaking terms are given by the usual collision 
integrals j( co ") in the equations for the diagonal components and by 

IV = 1.1 • 10~ 22 (y/x 5 ) MeV (42) 

for the non-diagonal component p^ = p*^. In the expression for T^ T we neglected the 
contribution from the majoron related processes. The effective potential V ab contains con- 
tributions from the usual weak interactions and from the neutrino-majoron interactions 
(26). We assume that the weak part is dominated by the charge asymmetric contribution, 
which is true in the case of a large charge asymmetry: 



1.5 Ci G F T 3 J dy ( Pab - p ab ) = B J dy (p ab - p ab ) . (43) 



with C\ defined in Eq. (6), dy = d 3 y/(27r) 3 , p is the antineutrino density matrix, and the 
coefficient 1.5 comes from (11/4) • 2((3)/it 2 , i.e. from normalization to photon number 
density. 

Equations (39-41) can be solved numerically but we can make reasonable estimates 
analytically in the following way (for more details see e.g. Refs. [1, 26]). In the case of 
strong coherence breaking, i.e. T^ T 3> H, Eq. (41) can be formally solved in the stationary 
point approximation i.e. putting r.h.s. equal zero. In this way we obtain: 

(/V ~ Ptt) _ P^ ~ Ptt („ gmfg W gV - C 2 Sm 2 32 /2E \ 

^ - - u TT - *r, T ~ v-^v7 T + ~ae 52 ) [ 1+ v„-v TT ) (44) 

where Sm 2 2 = m 2 — m 2 , s 2 = sin 28 and c 2 = cos 28. 

This expression could be substituted into Eqs. (39) and (40) for the diagonal components 
but one cannot obtain the closed equations for the latter because the potential V^ T contains 
an integral from p„ T over momentum, see Eq. (43). We assume for simplicity that only 

weak- interaction potential has noticeable off-diagonal components, i.e. V^ T ~ vffl . It can 
be true if e.g. |<7 rr | S> \g^\- To express V^ T through diagonal components we integrate 
Eq. (44) and similar equation for p^ T over momentum and subtract one from the other. To 
simplify the expressions let us present the diagonal part of the potential as a sum of charge 
symmetric (majoron) and antisymmetric (weak) parts: 

V m - V TT ee V™ + v[ w) (45) 

where the charge asymmetric part can be separated into two terms containing integrals 

from neutrino and antineutrino elements of the density matrix (see Eq. (43)), = V v —V v . 
After straightforward calculations we obtain: 



v^-i^-^r S2 s m 2 2 B 



4E 



Pw ~Ptt ( 1 ilV - 5ml 2 c 2 /2E 



vi x) + v v -v„\ V± x) + V U -V U 

Pnn ~ Ptt ( 1 if ixt - 6ml 2 c 2 /2E 



vi x) -v u + vA vi x) -v u + v v 



.(46) 
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This result is valid if |f| x) | > \V M \. In the case of the usual weak interactions considered 
e.g. in Refs. [7, 8] the situation is opposite, |V+ X ^| <C ivl^l, and the character of the 



approximate analytical solution is completely different. 
g l 

\Vfj, T \ < \s2dml 2 \/2E, as follows from Eq. (46). In this case we find 



For a rough estimate let us assume that |V+ X ^| > |V r i"'' > | and \V+^\ > \5m^ 2 \/2E, while 



V 4E J (y$ - v^y m 

The evolution of the difference of muonic and tauonic charges is determined by 

= f dy L<f) = j dy ( Pwt - p TT - p^ + p TT ) . (48) 

The collision integrals disappear from the time derivative of this difference since weak 
interactions conserve leptonic charges and majoron induced processes are assumed to be 
suppressed, see sec. 4 and in particular Eq. (32). The evolution of L is governed by the 
equation: 



dZ> T ) ( \g\ \" 4 / <5m§ 2 



dx V 1-5 x 10~ 7 J V 2.5 x 10~ 3 eV 2 




(49) 



and for the majoron coupling constant bounded by the conditions presented above, L^ T ) 
remains practically constant till nucleosynthesis epoch, i.e. x ~ 1. 

We have performed numerical calculations to solve Eqs. (39-41) in a similar way to 
Ref. [7]. The coupling constant matrix g a b is assumed to take the form: g a ^ = g5 aT 5t, T , for 
simplicity. The evolution of dimensionless chemical potentials and £ r for several values 
of \g\ are shown in Fig. 6. One can see that the lepton asymmetries of and v T equilibrate 
around T ~ 12MeV when \g\ =0. We have checked that the evolution remains the same for 
\g\ < 10~ 8 . As \g\ increases, the oscillations become less efficient and completely stop for 
\g\ ^ 10 -7 , which well agrees with the analytical estimates. It should be also noted that the 
lower bound on \g\ obtained from Eq. (49) does not depend on the magnitude of the initial 
lepton asymmetries. We have checked that the oscillations are similarly stuck for \g\ ^ 10~ 7 
even if the initial lepton asymmetries are taken to be smaller. The only difference is the 
temperature at which the asymmetries equilibrate when g = 0. 

We can repeat the same arguments for oscillations between v e and ty. The squared 
mass difference and mixing angle parameters for the oscillations between v e and v„i are 
5m,2i = m| — mf = 7.3 x 10~ 5 eV 2 and sin 2 9 = 0.315 [27]. The coupling constant matrix 
g a b is similarly approximated to be g a b = g 5 a ^i 5^ . Also the effective potential induced by 
the energy densities of electrons and positrons is taken into account, since here we consider 
oscillations including u e . The lower bound of \g\ becomes slightly relaxed due to the smaller 
mass difference, as can be seen from Fig. 7. The reason why the lepton asymmetries of u e 
and are not equilibrated completely when g = is that the mixing angle we adopt is 
not maximal. 

Lastly, we comment on the form of the coupling matrix, g^, for which the [y e — ty)- 
transformations are suppressed and, consequently, a large lepton asymmetry of the universe 
is allowed. So far we have assumed that g a b is diagonal in order to deal with both v T — 
and v e — zv oscillations similarly. However, in fact, a coupling matrix with more generic 
form works as well. We would like to stress that the equilibration of the muonic and tauonic 
lepton asymmetries is not harmful to our purpose, as long as these asymmetries are not 
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T [MeV] 



Figure 6: The evolutions of £a and £ T for several values of \g\ for the maximal mixing and 



bm\ 2 = 2.5 x 10- 



5 eV . The initial conditions are = —0.5 and ^ T = 1.0. 



erased through either oscillations or by the reactions shown in Figs. 3 and 4. Therefore 
we do not have to assume the hierarchical structure in g TT , g^r, and g^, and they are 
constrained only by the astrophysical bounds: 

| Sot | < 5 x 10~ 6 (50) 

with (a, b) = {(t, r), (//, r), (/i, //)}. In order to suppress the oscillations between v e and &y , 
at least one of these coupling constants should be larger than ~ 10~ 7 , while \g e e\i |<?eu|> 
and \g eT \ must be much smaller than ~ 10~ 7 . Thus we conclude that the coupling matrix 
suitable for our purpose should satisfy: 

-7 



10 



— 7 



< 



I See | j 

Max [ |(7 T1 



I fle/x | ; 1 9er \ 



< 10" 
<5 x 10" 



(51) 



In the simplest class of majoron models, the coupling matrix, g a b, is proportional to the 
neutrino mass matrix (m„) a f,. In particular, for the normal mass hierarchy, mi < mi < m^, 
the reconstructed mass matrix is often parametrized as [28] 

A 2 





( o 









g oc m u oc 





1 




.( 






1 







A 
A 




(52) 



with A ~ 0.2. Therefore the coupling matrix of this type can satisfy the conditions (51). Of 
course, in the more involved models, the coupling matrix is not necessarily proportional to 
the neutrino mass matrix. 



6 Discussions and Conclusions 

In the preceding sections we have shown that neutrino oscillations in primordial plasma 
can be blocked by an introduction of the neutrino-majoron interaction with a moderate 
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Figure 7: The evolutions of £ e and £„' for several values of \g\ with sin 2 9 = 0.315 and 
Sm^i = 7.3 x 10~ 5 eV 2 . The initial conditions are £ e = 0.1 and ^ = —0.5. 

value of coupling constants, |g| ^ 5 • 10 -6 . Here let us take up the case of \g\ ^ 10 -5 . For a 
definite discussion, we assume that the elements of coupling constant matrix whose indices 
involve electronic neutrino are much smaller than the other elements, i.e. \g e e\, \dea\ "C \g a b\ 
with (a, b) = {(r, r), (/x, r), (/i, //)}, so the electron type lepton number conserves during 
the relevant BBN epoch. As we have seen in section 4, the muonic and/or tauonic lepton 
number is not conserved for \g a b \ ~ 10 -5 . The fate of the large degeneracy in j/„ and v T 
depends on whether this L-violating interaction reached equilibrium before or after z/„ and 
v T decoupled from the thermal bath. 

First let us consider the former scenario. Since Vn and u T kept on to be in strong contact 
with the primeval plasma when the L-violating interactions were efficient, the muonic and/or 
tauonic lepton asymmetry vanished and thereby heated the plasma. At the same time 
majorons were abundantly produced and contributed to the effective number of neutrino 
species as AN U = 4/7. The CP asymmetric part of the effective potential of neutrino due 
to the muonic and tauonic lepton asymmetries vanished, but there would be a similar (but 
possibly smaller) term due to the charge asymmetry of v e (see Eqs. (7) and (8)). Also there 
exists the contribution the the effective potential due to abundant majorons in the plasma 
(see Eq. (25)). A further important point is that L-violating processes, shown in Fig. 3, 
induces the effective potential between neutrino and antineutrino, since (cvb) no longer 
vanishes. However, neutrino oscillations between v e and can be treated in the same way 
as before, since electron type lepton number remains a well conserved quantity. Therefore, 
we can deduce from the previous results that neutrino oscillations between v e and are 
blocked and £ e remains intact. It should be noted that equilibrium majorons, instead of 
excessive Un and u T , contribute to the extra radiation, which speeds up the expansion rate, 
and that its contribution, AN U = 4/7, well satisfies the constraint shown in Eq. (3). 

Next we discuss the other possibility that the L-violating interactions come to thermal 
equilibrium after decoupling of and v T . Majorons are also produced, as in the previous 
case. However, since and v T cannot exchange energy with the primeval plasma, the 
total energy of majorons, Vn and v T must be conserved. Therefore the excessive energy 
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previously stored in i/„ and v T is redistributed among majorons, u^, and v T , and their 
distribution would be like Eq. (33) but without any charge asymmetry, i.e. with £ = £. It is 
thus clear that AiV^ remains unchanged before and after L- violating interaction comes into 
equilibrium. Our concern here is whether the modified distribution of ty affects that of v e 
through neutrino oscillations. However, the oscillations are blocked on the same ground. 

In connection with the second case, there is an interesting possibility that the majoron- 
exchange reactions: v^ T v^ T v e u e , v^ T v^ T <-»■ v e v e , v^ T v^ T <-> P e P e , v^v^ <-»■ v e v e and 
^h,t^u,,t & e P e , could produce additional electronic neutrinos and antineutrinos. These 
processes can proceed if \g ee \ and \g ea \ with a = [i, t are sizable ^ 10~ 5 . If this is the 
case, the relative abundance of v e and P e with respect to photons and electrons would 
be larger than in the standard model. This would lead to a later neutron-proton freezing 
and to a smaller number density of survived neutrons, which might compensate the effect of 
additional energy density stored in all species of neutrinos and majorons and even overshoot 
it. Then we do not need to suppress neutrino oscillations, therefore the coupling constant 
matrix is allowed to take rather arbitrary values as long as the L-violating processes 
come into equilibrium after decoupling of muonic and tauonic neutrinos. 

Thus our examination indicates that a conspiracy between the speed-up effect and the 
shift of (5 equilibrium is possible for \g a b\ > 10 -5 , as long as \g ee \, |<7ea| ^ |<7a&| with a, b = 
H, t. Furthermore, even for \g ee \, \g e a \ ~ 10~ 5 , such a cancellation might be possible in a 
somewhat different way. This considerably extends the allowed parameter space, although 
rigorous study might be necessary to obtain further quantitative results. 

In this paper we have shown that the hypothetical neutrino- majoron interaction can 
suppress neutrino oscillations in the primordial plasma to prevent lepton asymmetries of all 
neutrino species from being equilibrated. The exact form of the effective potential induced 
by this interaction is calculated. We have found an allowed range of the coupling constant: 
10~ 7 < \g\ < 5- 10~ 6 (more precisely, see Eq. (51)), which satisfies the astrophysical bounds 
and makes the scenario operative. For the coupling constant in this range, u e — u^i oscillation 
in the early Universe is blocked, thereby keeping the cosmological lepton asymmetry of 
electron type unchanged. The upper bound comes from the requirement that lepton number 
is effectively conserved, and the lower bound is obtained from the study of the evolution of 
the lepton asymmetries both analytically and numerically, in two flavor approximation. The 
constant matrix in the simplest class of majoron models can satisfy the desired constraints, 
in the case of the normal mass hierarchy. Furthermore, even for \g\ ^ 10~ 5 , the large energy 
density necessary to cancel the effect of £ e can be supplied by majorons, v„ and v T , even if 
muonic and/or tauonic lepton asymmetries were erased. Thus we conclude that an addition 
of the majoron field to the standard model can reopen a possibility that the effect of £ e 
is compensated by large £^ iT (or by the extra energy of majoron itself), thereby curing a 
probable discrepancy between the BBN and CMBR. 

A Derivation of effective potential due to majoron-neutrino 
interaction 

Here we present the derivation of the effective potentials Eqs. (23) and (25). First the nota- 
tions and conventions that we adopt here, are listed. The metric is taken to be (— , +, +, +). 
The gamma matrices in the chiral representation are 

'--'(?;)• '~'U o)< <53> 
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where Cj are the Pauli matrices. The gamma matrices in the Dirac representation arc 
related to those in the chiral representation as 



7b: 



Dirac — ^7chiral^ ' ^ 



1/1 1 



V2 \ i -i r 



(54) 



and they are given by 



7 



1 



, 7 = i 



a; 



v -1 J ' ' ^ -<Ti 

The charge conjugation matrix C is defined as 

_.( (7 2 

-(o 

Y ^2 \ 



75 



1 

1 



(55) 



(chiral) 



(Dirac) 



(56) 



Hereafter we work with the gamma matrices in the Dirac representation. The momentum 
expansion of a free Dirac field ip(x) with mass m is given by 

^(x) = I d P ( a (P' S MP' s)e- iE * t+ip * + &t( p> a ) u ( Pj s ) e iE p*- i P x ) , 

s J 

= E / d P («(P. S MP' s ) + bf ("P. S M"P' s )) e " P ° i+ ' PX . ( 57 ) 



where we defined p° = epE'p = e p -\/p 2 + "i 2 with e p = ±1. The sign of e v is chosen so 
that the positive and negative-energy solutions are reproduced, s = ± describes the sign of 
the spin projection eigenvalue, and we have used the notation dp = d 3 p/(27r) 3 . The Dirac 
spinors u(p, s) and v(p, s) take the form, 



u(p,s) 



v(p,s) = 



E p + m 



^2E p (E p + m) V ~ p 



-pa 



2E p (E p + m) V E P + m 



Xu(s), 



Xv(s)- 



(58) 



The two-component spinors Xu(s) and Xv(s) are not independent but related to each other 
through: Xv(s) = — i&2Xu(s)*- If we choose them to be eigenstates of the helicity operator 
p • ct/IpI) they are 



, / e"^ 2 cos § \ 

X«(+) = \ P if/2 sin i ) . XuH = 



_ e ~i<p/2 gin | 

e^cosf 



X«(+) = 



-^/ 2 sin| 



D iip/2 



cos 



. Xv(-) 



-if/ 2 



COS 



_e<v/2 S i n : 



(59) 



where 6 and ip are defined as p = |p| (sin cost/?, sin sin </?, cos 0). For flipped momentum 
— p, {6, ip} should be replaced with {ir — 6, ip + ir}. 
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In deriving the effective potential, the neutrino field u(x) is approximated to be a mass- 
less left-handed field : 

Vi (x) = j dp(a i (p) Up + b l (-p^v_ p )e- i P 0t+l ^, (60) 



where 



1 + 75 , v i / n , . 

Up = — u(p,-) = -j= I 1 I xu(-), 

«- P = — 2~ «(-P>-) = ^ I j j Xu(+), 
1 + 75 



2 



-Kp>+) = T^f ( 1 ) x "(")' 



f-P ^ l Ap L v{-p,+) = ^=[ \ j Xu (+). (61) 

The annihilation and creation operators a, 6, a) and 6^ satisfy the anticommutation relations. 

{(Ufa), a,-(q)t} = ^^-^(p-q), 

{6i(p), 6,-(q)t} = ^^-^(p-q), 

others = 0. (62) 

Similarly, the momentum expansion of the free majoron field is 

X (x) = jdp(a x (p)u x (p)+a x (-p)h x (-p))e^° t+ ^, (63) 

where 

%(P) = «x(P) = "7^" ( 64 ) 
The annihilation and creation operators a x , a x satisfy the commutation relations: 

a x (p), a x (q)t] = (2vr) 3 ^ (p - q), 

others = 0. (65) 

Next we derive the effective Hamiltonian Eq. (20). The equation of motion for the 
majoron field \ is 

d^ X = ~\ (gab v T a Cv h + g* ab v\Cvl) . (66) 
To solve this equation, we define the majoron and neutrino fields in momentum space : 



d 4 p 



*a(x) = J ^ji? a (p) e **. (67) 



(2VT) 
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Then the solution of Eq. (66) is 

X(P) = ^{^{gabUp-qfCh^+glbhi-p + q^CDli-q)}. (68) 



Thus the effective Hamiltonian is: 



-eff ' 



d 4 P 



d A q 

d 4 r 



(2vr) 4 2p 2 
+g: b h{-p + q) ] CDl{-q)} J 
x {g cd D c {r) T Cv d {s) + g* cd D d {-s) ] CD*{-r)) . 



(27T) 



{gabMp- q) T Ch(q) 

d As isx 
(2vr) 4 



(69) 



where we substituted Eq. (68) into the second equation. Note that the numerical factor 
1/2 is inserted in front of the last equation. In a first-order perturbative approximation the 
neutrino and majoron fields can be set to be free fields: 



h(p) 
x(p) 



2ir6(p° -e p E p )vi(p), 
2tt5( P ° -e p E p ) X (p), 



(70) 



where 



X(P) = o x (p)u x (p) + a x (-p) t t> x (-p). 



(71) 



Thus the effective Hamiltonian becomes 



-J dpdqdrds (2^) 3 ^(p + q + r + s)e^ tot * 

F(p,q)F(r,s) 
e p e q E v Eq(\ - e p e q cos pq ) ' 



(72) 



where 



Etot 
cos 6> pq 



p q 



IpI |q| 



F(p,q) = g a b v a (p) T C v b (ci) + SW ) ffc(-q) t CV*(-p). 



(73) 



Once the effective Hamiltonian is obtained, we can calculate the effective potential for 
neutrinos according to Ref. [23] . All we have to do is to evaluate 



Weff(t = 0),a,-(p)ta i (p)l) = -(2vr) 3 5( 3 )(0)[y p , / 9 p ] ij , 



(74) 



where [Vp] ab is the effective potential matrix. Using the anticommutation relations (62), we 
obtain the final result: 



ab 



-I 



dq 



1 



4 IpI |q| 



ab 



(75) 



19 



Next we derive the effective Hamiltonian, which describes neutrino- majoron scattering. 
In this case we integrate out the neutrino field instead of majoron field. The equation of 
motion for the neutrino is 



-nV^a + *x<?:&cv 6 * = o. 

d 4 q 



The solution is then 

*«(P) = |^<4 J (0X(P + «)7"7 O ^( ? ). 
Substituting this result into Eq. (19), we obtain the effective Hamiltonian: 



(76) 
(77) 



n (ux) (t) - - l -a,at f d 3 x f ^ ^ ^ c i(p+(,+r+s)x 
n eS \t) - 2 g ab g bc j a xy ( ^ )4 ( ^ )4 ( ^ )4 ( ^ )4 e 

xx(p)x(?)^(-r)y5.(s) (</ + '" ] " 



ab 



(q + r) 2 ' 

/" tip dq dr ds (2tt) 3 5 (3) (p + q + r + s) e - iEtott 



xx(p)x(q)^a(-r)7 M ^(s) 



2g • r 



g =egS q ,rO=e r £; r 



(78) 



where we take both neutrino and majoron fields to be on-shell in the second equation. 
Applying Eq. (74), we obtain the effective potential due to the neutrino- majoron scattering: 

/x(q) 



ah 



dq 



4|p| |q| 



9 ] 9 



ab 



where the number density of majorons with momentum p, is defined as 



4(p) a x (p' 



(79) 



(80) 



B Invariant amplitude squared for 2v <-> 2v 

Here we calculate the invariant amplitude squared for the reaction, 2v <-> 2v. Provided that 
the diagonal coupling constants are larger than the non-diagonal ones, and that one of the 
diagonal coupling constants dominates, we consider neutrinos of the flavor with the largest 
coupling. Hereafter we drop the sub-indices for flavor. The diagrams, which contribute to 
the reaction, are shown in Fig. 3. 

The S-matrix and the invariant amplitude squared are related as 

Sp a = 27r(-i)5 {4) {p l3 -p a )Mf3a, 
|M inv | 2 = (27r)- (i l[[(27Tf2E a \l[[(27Tf2E f3 \\M f 3 a \ 2 , (81) 

where a and (3 represent symbolically the initial and final states . Each matrix element for 
s-, t- and u-channel diagrams is given as 

q 2 1 

M pLpiP2 = (2?r) 3 {pi +p2) 2 ( U li Cu Pl)( V P4 Cv P3)> 

t 9 2 1 T T 

"^P3P4,P1P2 = _ (27r)3 ( pi _ p3 )2^P3 C ' U Pl)( V P4 C '' U P2)' 

M1, P1P2 = j£y 3 {pi l p4) 2 (< Cu M 3 Cu P2 ), (82) 
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where the four momenta satisfy the conservation condition, p\ + p 2 = p 3 +P4- If we use the 
explicit expressions for u p and v p derived in the previous section, we obtain 

„2 i o 

A/f( s ) 4- A/f(*) A- M^> = — — - (83^ 

J l P3P4,PlP2^ J l P3P4,PlP2^ J l P3P4,PlP2 (27r) 3 ^E 1 E 2 E 3 E 4: 4' 

Substituting this result into Eq. (81), the invariant amplitude squared becomes 

|M inv | 2 = 9|<?| 4 . (84) 

C Formulas for the collisional integral in the Boltzmann ap- 
proximation 

Here we write down the formulas for the collisional integral in the Boltzmann approximation 
for completeness. For the derivation, see Ref. [29]. 

-e- 2 *-^ 4 ) (yi + -2/4) = 47r 2 e-f\ 



rd 3 y 2 


fd 3 y 3 f 


1 V2 J 


V3 J 



r in. ( *zi J r fi± e -„- msW{m + m _ _ , _ 32lr3 . (85) 

J Vl J V2 J V3 J V4 
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